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quantities in the region of the Large Mixing Angle-MSW solar solution - order of 
(— I \ magnitudewise by only six adjustable parameters (Higgs vacuum expectation values) 

to also give an agreeing prediction for the amount of baryogenesis produced in the 
early time of cosmology. We use Fukugita-Yanagida scheme and take into account 
now also the effect of the remormalisation equation for the Dirac neutrino sector 
' from Planck scale to the see-saw scale. The present version of our model with many 

approximately conserved (gauge) charges distinguishing various left- and right- Weyl 
particles has the largest matrix element of the mass matrix for the three flavour 
see-saw neutrinos being the off-diagonal elements associated with the second- and 
third-proto-flavour and gives the ratio of baryon number density to the entropy 
density to 2.59+2^25 ^ ^O"" which agrees perfectly well as do also all the fermion 
masses and their mixing angles order of magnitudewise. 
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1 Introduction 



Recently we have improved a rather specific model seeking to fit/explain quark and 
lepton masses and mixing angles using as the reason for the large mass ratios and usually 
small mixing angles, except two neutrino mixing angles, approximately conserved quan- 
tum numbers invented for that purpose into our model. The approximately conserved 
quantum numbers - which are really assumed to be gauged, although that may not be 
so crucial - are supposed to be broken by various Higgs fields. We could get a good fit 
of all the masses and mixing angles. However, this previous version of our model fails 
by providing too little baryogenesis to be produced in the early Big Bang. To calculate 
the baryogenesis we use the Fukugita-Yanagida scheme for producing a. B — L excess at 
the time when the temperature of the Universe passes the scale of the see-saw neutrino 
masses; in fact the see-saw mechanism for neutrinos masses is incorporated into our model 
in as far as we have three right-handed neutrinos having Major ana masses acquired via a 
certain Higgs vacuum expectation value (VEV), {4>b-l)- 

In such model(s) we can play around with the system of charges, but most importantly 
we have played with developing the charge assignments of the Higgs fields breaking the 
(gauge) group so that these charges are only approximately conserved. It is the detailed 
quantum numbers of the speculated Higgses that we used to develop to fit more and 
more experimental informations. In the last improvement we got the model to predict 
Large Mixing Angle-MSW (LMA-MSW) solution rather than Small Mixing Angle- 
MSW (SMA-MSW) solution by replacing the two of our Higgs fields which are able to 
cause transitions between first and second generation quarks and leptons by a new pair 
which actually turned out to have more elegant quantum numbers. 

As we are now fitting all the masses and mixing angles, the constraints on our model 
are so tight that we hardly can change it anymore. However, we shall see below - and 
that is the main point of the present article - that we could still change the quantum 
numbers of the Higgs field delivering the over all scale of the see-saw neutrino masses. 

The difficulty of getting enough baryons with our type of model hangs together with 
that even the largest element in the (Dirac) neutrino matrix is suppressed compared to 
unity by a factor 10~^. Only our up-type quark and right-handed neutrino mass matrix 
have order of unity elements. Models with two Weinberg-Salam Higgs doublets can on 
the contrary have a tan/3 that is large and give suppressions to the down-type quark 
and charged lepton masses so that the Yukawa couplings for these (lepton and down-type 
quark) particles could be as large as of order unity. That is to say two Higgs doublet 
models may in general tend to get stronger Yukawa couplings and thereby easier CP 
violation since the latter comes only in one-loop accuracy. 

It means that we are in need for possible enhancements for baryogenesis. The possi- 
bility for that which it actually turns out that we can obtain very elegantly in our model 
is to make the mass difference between some of the see-saw neutrinos smaller than in the 
previous model. In fact the CP violating asymmetry parameter is essentially inversely 
proportional to the mass difference between the see-saw neutrino producing the asymme- 
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try in its decay and another see-saw neutrino. Due to the self-energy contribution almost 
degenerate neutrinos can even enhance the asymmetry by orders of magnitude compared 
to what it would be with the right-handed neutrino mass ratios of order unity. It would 
therefore be a progress in the direction of producing enough baryons to put two of our 
see-saw neutrinos almost degenerate in mass. It is very natural to obtain this situation 
by letting an off-diagonal element in the see-saw neutrino mass matrix dominate. That 
may be done by adjusting the quantum numbers of that Higgs field (f)B-L which gives the 
masses and the mass scale of the see-saw neutrinos so that it causes the transition from 
one flavour to another. 

This article is organised as follows: in the next section, we present general assumptions 
for the model building, and we review our model - the family replicated gauge group 
model. In section 3 we define our notation for the fermion mass terms in the Lagrangian 
and their mixing angles. Then, in section 4 the renormalisation group equations of all 
sectors are presented. The calculation is described in section 5 and the results for the 
fermion masses and mixing angles are presented in section 6. In section 7 we discuss the 
Fukugita-Yanagida scheme of baryon number production, and in section 8 the wash-out 
of the produced B — L excess is approximated and our results of the calculation are given 
in section 9. A discussion on the proton decay goes into section 10. Finally, section 11 
contains our conclusions. 

2 Model with many mass protecting charges 
2.1 General feature 

The main point of the model we use here is that the "small hierarchies" of the quark 
and lepton masses are due to approximately conserved quantum numbers [§] preventing 
the masses from being different from zero, i.e., that the Yukawa couplings which are 
in Standard Model mysteriously small should really be understood as being so due to 
that they need breaking of some charge conservations. In fact we consider these Yukawa 
couplings just some effective couplings representing at the fundamental level more com- 
plicated vertex diagrams involving attachments to the VEVs of Higgs fields breaking the 
charge conservations in question (see Fig. |1|). 

We summarise here our philosophy to build a model that can predict, rather say fits 
by it having the all fermion quantities including not only the neutrino oscillations but 
also baryogenesis in following ingredients: 

• All fundamental coupling constants and masses are of order unity. 

All gauge and Yukawa couplings are of order one at the Planck scale, in fact, we 
take the fundamental scale as the Planck one. From this point of view we treat 
them as random numbers, and partly to be able to "explain" CP violation 
we assume that all Yukawa couplings are also complex random numbers. Also 
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the masses of particles should be of order Planck scale unless they are mass 
protected, so as to have masses determined by the VEVs of Higgs fields say. 

• Only VEVs are allowed far from order one. 

Thus also the scalar field masses corresponding to the VEVs are allowed to 
be very small, since otherwise it would be inconsistent to have = A {(pws) 
unless A has huge value. 

This assumption could gain supports from 

(a) the well-known weak scale being extremely small compare to the Planck scale. 

(h) dynamical symmetry e.g. Nambu-Jona-Lasinio can lead to equation for the 
logarithm of the VEV scale. 

(c) ui superconductivity VEVs tend very small compared to the atomic dimension- 
ality argument. 

• There are a priori several conserved charges and the various Weyl components of 
quarks and leptons have usually different charges, i.e., quantum numbers. 

Since we let Weyl components of different families have in general different 
quantum numbers of these kinds, they are to be considered horizontal, but 
they should better not bring together into multiplets different families because 
that would make the remarkably large inter-family mass ratios difficult to in- 
corporate (unless the group is strongly broken). 

• Incorporate the see-saw mechanism |^. 

A new scale must be introduced into the model. In our case we assume a new 
scale about 10^° GeV. This scale becomes crudely the masses of right-handed 
(Majorana) neutrinos, and therefore B — L quantum charge being broken at 
that energy scale. 




Figure 1: Diagram giving the effective Standard Model Dirac Yukawa coupling. Here 
Weinberg-Salam Higgs field is denoted as (pws- 
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Note that we do not assume the existence of Supersymmetry or Grand Unification. How- 
ever, our prediction came out that SU{5) mass relations are order of magnitudewise vahd 
due to the diagonal elements of different mass matrices need same quantum numbers 
being broken in our specific model. 

2.2 Our specific model of many conserved charges 

We have already investigated a model || which can predict all quark and charged 
lepton quantities including baryogenesisQ This model is of the type described in the 
foregoing subsection with the further good restriction that the quantum numbers for the 
Weyl particles are so chosen that they obey the anomaly cancelation conditions required 
for them to be imagined gauged charges. It is namely easily seen that our gauge group 
has no anomalies - neither gauge nor mixed ones - because it has the Standard Model 
group plus a. B — L separately for each family. In fact our model has the specific group 



where SMGi = SU{3)i x SU{2)i x U{l)i denotes the Standard Model gauge group for 
each generation {i = 1, 2, 3); x means the Cartesian product. However, it is unlikely to be 
important to have this special group, provided we have sufficiently many conserved charges 
separating in quantum number the various Weyl fields. Actually the important thing is 
that one by means of the quantum numbers can construct some Higgs fields as we propose 
them for the breaking of the group to the Standard Model with much the same powers 
of fields being needed in the various mass matrix elements. The latter may in fact not be 
difficult to obtain with other groups, at least the non-abelian part of our proposed group 
is totally irrelevant for the fitting. 

We should emphasize that due to the non-zero neutrino masses using the see-saw 
picture, it is necessary to introduce a right-handed neutrino, or preferably several - we 
have three in our model - having Majorana masses. In our model we associate the see- 
saw mass scale with a gauged B — L charge. We let the latter come from the diagonal 
subgroup of the Cartan product of a f/(l)s_i„, {i = 1, 2, 3) for each family. 

Note that this family replicated gauge group, eq. (|ID, is the maximal gauge group 
under the following assumptions: 

(1) Considering only part of the gauge group acting non-trivially on the known 45 Weyl 
fermions of the Standard Model and the additional three heavy see-saw (right- 
handed) neutrinos, i.e., our gauge group is a subgroup of f/(48). 

(2) Avoid gauge transformations transforming a Weyl state from one irreducible rep- 
resentation of the Standard Model group into another irreducible representation, 

^The model provided only the SMA-MSW solution. A main point of present article is to calculate 
baryogenesis with LMA-MSW solution predicting mass matrices including running effects of all sectors 





(see sections 0, 7-9). 
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i.e., there is no place for the quarks and leptons occupying the same irreducible 
representation under our gauge group. 

(3) No gauge nor mixed anomalies due to the renormalisability requirement. 

2.3 Description of quantum numbers of mass protected particles 
in model 

In addition to the practically any non-mass-protected particles, which we have already 
assumed to be present with Planck scale masses, we have as is easily seen some mass 
protected particles. The Weyl particles without chiral partners in our model are the well- 
known Weyl components of the quarks and leptons in the Standard Model plus three see- 
saw (right-handed) heavy neutrinos, which are all easily seen to be mass protected by our 
gauge group eq.(|I|). These 48 Weyl particles have lower mass than the Planck scale. They 
fall into three families each consisting of the 16 Weyl particles of a usual Standard Model 
generation plus one see-saw particle. In this way we can label these particle as proto-left- 
handed or proto-right-handed w-quark, (i-quark, electron etc. To get the quantum numbers 
under our model gauge group for a given proto-irreducible representation, we proceed in 
the following way: We note the generation number i {i = 1, 2, 3) of the particle for which 
we want quantum numbers and we look up, in the Standard Model, what are the quantum 
numbers of the irreducible representation in question and what is the B — L quantum 
number. Then we take trivial quantum numbers for the two generations j different from 
i and put the index i to the gauge group's names, so we put the representations of the 
groups S'f/(3)j, SU{2)i and U{l)i equal to the just asked quantum numbers. 

For instance, if we want to find the quantum numbers of the proto-left-handed bot- 
tom quark, we note that the quantum numbers of the left-handed bottom quark in the 
Standard Model are weak hypercharge y/2 = 1/6, doublet under SU{2) and triplet under 
SU{3), while B — L is equal to the baryon number = 1/3. Moreover, ignoring mixing 
angles, the generation is denoted as number i = 3. The latter fact means that all the 
quantum numbers for SMGj j = 1,2 are trivial. Also the baryon number minus lepton 
number for the proto-generation number one and two are zero: only the quantum numbers 
associated with proto-generation three are non-trivial. Thus, in our model, the quantum 
numbers of the proto-left-handed bottom quark are 1/3/2 = 1/6, doublet under SU{2)3, 
triplet under SU^S)^ and {B — L)^ = 1/3. For each proto-generation the following charge 
quantisation rule applies 

1 + 1 + 1 = (modi) , (2) 

where ti and di are the triality and duality for the z'th proto-generation gauge groups 
SU{3)i and SU{2)i respectively. 

Combining eq. (^ with the principle of taking the smallest possible representation of 
the groups S'f/(3)i and SU{2)i, it is sufficient to specify the six Abelian quantum numbers 
yi/2 and {B — L)i in order to completely specify the gauge quantum numbers of the fields, 
i.e. of the Higgs fields and fermion fields. Using this rule we easily specify the fermion 
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representations as in Table p.3| . However, as already mentioned in section 2.2 these non- 
abelian SU{3)i and SU{2)i are really not important for our fit since they just follow the 
weak hypercharges like "slaves". 

Note that fermion quantum numbers for each proto-generation gauge group SMGi x 
U{X)B-L,i are obtainable as considering it a subgroup of 5*0(10), i.e. our gauge group 
eq. (P is really a subgroup of S'O(IO)^. However, we avoid the not well-agreeing accurate 
predictions of simplest grand unification by not having neither SU{5) 's nor SO{10) 's: the 
mass spectra, and also proton decay. 
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Table 1: All U{1) quantum charges for the proto-fermions and Higgs fields in the model. 
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2.4 Higgses breaking the family replicated gauge group to the 
Standard Model 



The gauge group ^^^2,3 {SMGi x f/(l)i3_L,i) is at first spontaneously broken down at one or 
two orders of magnitude below the Planck scale, by 5 different Higgs fields, to the gauge 
group SMG X U{1)b~l which is the diagonal subgroup of the original one: 



X 

1=1,2,3 



{SMG, X f/(l),_,„:) -> ( J,, {SMGi X , = SMG x U{1)b-l ■ 

V ' / diagonal 

(3) 

This diagonal subgroup is further broken down by yet two more Higgs fields — the 
Weinberg-Salam Higgs field (pws and another Higgs field — to SU{3) x U{l)em- 

The VEV of the 4>b-l Higgs field is taken to be about 10^° GeV and is designed to break 
the gauged B — L quantum number. In other words the VEV, gives the see-saw 

scale (All U{1) quantum charged for the proto-fermions and Higgs fields are presented in 
Table 



We should emphanise here that the quantum numbers of the (ps-L has been changed 
relative to our previous publication due to the baryogenesis point of view. The point is 
that baryogenesis calculated the way described below with the "old" quantum charge of 

4>B-L- 

Q^B-L I ^^^ = (0,0, 0,0, 0,2) (4) 

would produce insufficient amount of baryon minus lepton number, whereas we shall see 
that quantum numbers making the off-diagonal elements (2, 3) and (3, 2) of the right- 
handed neutrino mass matrix be the dominant ones as in Table |2l 



(0,0,0,0,1,1) . (5) 



3 Fermion masses and their mixing angles 

In this section we define our notations for the mass terms of the charged fermion sectors 
and the neutrino sectors. These define also the mixing angle unitary matrices for quark 
sector and neutrino sector, respectively. 



3.1 Charged fermion masses and their mixing angles 

The full Lagrangian of our model is that of a gauge theory with scalars and Weyl 
particles of which many are the ones at the Planck scale with which we do not go into 
details. We shall therefore leave it to the readers imagination and only present the notation 
for the Weinberg-Salam Yukawa part effective Lagrangian of the Standard Model form 

-^charged-fermion— mass ~ Q l^u^WsUr ~\~ Ql^D^wS^R ~^ -^L^B*^ ws-^'fl "I" h.C. (6) 
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Here ^ws is the Weinberg-Salam Higgs field, Ql denotes the three SU{2) doublets of 
left-handed quarks, Ur denotes the three singlets of right-handed up-type quarks and 
Yu is the three-by-three Yukawa coupling matrix for the up-type quarks. Similarly Yd 
and Ye are the Yukawa coupling matrices for the down-type quarks and charged leptons 
respectively. The SU (2) doublets ^ws and Ql can be represented as 2 component column 
vectors and we then define: 



and 



where are the CP conjugates of the three left-handed up-type quarks. 

The effective Yukawa couplings Ye, for example, are obtained from identifying the 
lowest order vertex in eq. with a vertex Fig. |I|. 

After electroweak symmetry breaking the Weinberg-Salam Higgs field gets a VEV and 
we obtain the following mass terms in the Lagrangian: 

-^charged— fermion— mass 

Ur^ Mu Ua + Mr, Dr, + Er Me E« + h.c. (9) 

where the mass matrices are related to the Yukawa coupling matrices and Weinberg-Salam 
Higgs VEV by: 

M = r%^. (10) 

v2 

We have chosen the normalisation from the Fermi coupling constant: 

= 246 GeV . (11) 

In order to obtain the masses from the mass matrices, M^, and M^, we must diago- 
nalise them to find their eigenvalues. In particular we can find unitary matrices, for 
the up-type quarks, for the down-type quarks and Ve for the charged leptons: 

V^M^MlVu = diag (m2,m2,m?) , (12) 

V^MnMlVn = diag {ml ml ml) , (13) 

V^MeMIVe = die.g [ml, ml, ml) . (14) 

The quark mixing matrix is then defined with these unitary matrices as 0: 

VaKM = V^Vn . (15) 

This VcKM unitary matrix is parameterised as follows: 

/ C12C13 S12C13 si3e""^" 

VcKM = ~'Sl2C23 " Ci2S23Si3e*'^^^ C12C23 — Si2S23Si3e*'^^^ S23C13 

\ ■Sl2'S23 ~ Ci2C23Si3e*'^" " C12S23 — 5120235136**^^^ C23C13 

/ Vud Vus 

Vcs Kb I , (16) 



V Vt 



ts 
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where Cij = cos% and Sij = smOij for the generation labels z, j = 1, 2, 3 and 5i3 is a CP 
violating phase in the quark sector. 



3.2 Neutrino masses and their mixing angles 

The left-handed neutrinos are massless in Standard Model because of non-existence of 
right-handed neutrinos and the Weinberg-Salam Higgs field having only weakhyper charge 
y/2 = 1/2, and not |//2 = 1 as is required to give Majorana neutrino a mass. Thus, 
in order to let left-handed neutrinos be massive particles as is very strongly indicated 
by experiments, we assume the existence of the three very heavy right-handed neutrinos, 
which are mass-protected so as finally get masses from the "new" Higgs field, (pB-L, at 
an energy scale of about 10^° GeV in our present model. 

We use the gauged B — L charge to mass-protect the right-handed neutrinos; in fact 
we use the total - diagonal - one in fact we break U{1)b-l,i x U{1)b-l,2 x f^(l)_B-L,3 
D U{1)b~-l at a much higher energy scale, say about 10^^ GeV. 

The assumption of the existence of three right-handed Majorana neutrinos at a high 
scale gives rise to the addition of Majorana mass terms to the Lagrangian: 

= ]^{nzyMnL + h.c. (17) 



where 




Mr M.^ 



Here is the left-right transition mass term - Dirac neutrino mass term - and and 
Mr are the isosinglet Majorana mass terms of left-handed and right-handed neutrinos, 
respectively. 

Due to mass-protection by the Standard Model gauge symmetry, the left-handed Ma- 
jorana mass terms, M^, are negligible in our model with a fundamental scale set by the 
Planck mass. Then, naturally, the effective light (left-left transition) neutrino mass matrix 
can be obtained via the see-saw mechanism lEi: 



^eff ■ 



■.M^M],\M!^f . (19) 



In the framework of the three active neutrinos model, the flavour eigenstates {a = 
e, z/, r) are related to the mass eigenstates z/j {i = 1,2,3) in the vacuum by a unitary 
matrix Vmns, 

Wa) = Xl(^MNs)m Wi) ■ (20) 
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Here Kins is the three-by-three Maki-Nakagawa-Sakata (MNS) mixing matrix 
is parameterised by 



( C12C13 

— 5i2C23 — Ci2S23'5l3e_ 
V 512323 - Ci2C235l3e*'''^ 

/ e*^ 
X e*'^ 
V 1 



S12C13 

C12C23 512523^136 



i(5i3 



-C12S23 - Si2C23Si3e 



si3e 
523C13 
C23C13 



which 



(21) 



where Cij = cos 9ij and Sij = sin 9ij are the neutrino mixing parameters and ^13 is a 
neutrino CP violating phase. Note that, due to the existence of Majorana neutrinos, we 
have two additional CP violating Majorana phases ip and ip, which are also included in 
the MNS unitary mixing matrix. 

In order to get predictions for the neutrino masses from the effective mass matrix, 
Mefj, we have to diagonalise this matrix using a unitary matrix, V^s, to find the mass 
eigenvalues: 

KffMeffMjgVjff = diag(m?, mj, ml) . (22) 

With the charged lepton unitary matrix Ve, eq. (|T^) , we can then find the neutrino mixing 
matrix: 



Ve Kir 



(23) 



Obviously, we should compare these theoretical predictions with experimentally measured 
quantities, therefore we define: 





= - 


- ml 




= ^l- 


-ml 


tan^ 9q 


= tan^ 


9l2 


tan'^ ^atm 


= tan^ 


923 


tan^ ^chooz 


= tan^ 


9l3 



1 5 



(24) 
(25) 
(26) 
(27) 
(28) 



Note that since we use the philosophy of order of magnitudewise predictions (see 
section ^ with complex order one coupling constants, our model is capable of making 
predictions for these three phases, the CP violating phase 613 and the two Majorana 
phases; put simply, we assume that all these phases are of order 7r/2, i.e. essentially 
maximal CP violations. 



4 Renormalisation group equations 

From the Planck scale down to the see-saw scale or rather from where our gauge group is 
broken down to SMC x U{1)b~l we use the one-loop renormalisation group running of 
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the Yukawa coupling constant matrices and the gauge couphngs including the running 
of Dirac neutrino Yukawa coupling: 



^dg2 19 3 



ISvr^^ = -7^73^ , (31) 

167r^^ = ^ {Yu{Y,y - r,(r,)t) + { n - (11^? + ^gl + 8gl) } , (32) 

levr^^ = I {YniY.y - YAYu)') Y, + { - (i,? + + 8,3^)} Y, , (33) 

le^r^^ = I (r.(r.)t - y;(y;)t) + 1 _ (^^^ _^ 9 ,x | ^^^^ 



IGvr 



^ = ^ (k(K)^ - r.(F.)t) K + { - (^,? + ^,,^)} K , (35) 

n = Tr(3y;tr^ + 3y;ty^ + yjy^ + ytyj , (36) 



where t = In and /i is the renormalisation point. 

However, below the see-saw scale the right-handed neutrino are no more relevant and 
the Dirac neutrino terms in the renormalisation group equations should be removed, and 
also the Dirac neutrino Yukawa couplings themselves are not accessible anymore. That 
means that, from the see-saw scale down to the experimental scale (1 GeV), the only 
leptonic Yukawa /3-functions should be changed as follows: 

16vr^^ = \ {Y^iY,)') + { Fs - (^^? + {gl) ] Y, . (37) 

Note that the quantity, Yg, must be also changed below the see-saw scale: 

y5 = Tr(3FjF^ + 3rjF^ + FiF^) . (38) 



Further, we should evolve the effective neutrino mass matrix considered as a whole 



running as an irrelevant or nonrenormalisable 5 dimensional term [|I0| from the see-saw 
scale, set by {(pB-i) to our model, to 1 GeV: 



167r' 



,dM, 



eff 



dt 



[-3gl + 2X + 2Ys) - - {m,s {Y,Y^) + {Y.Y^f M,^) , 



(39) 



where Yg defined in eq. ( |38| ) and in this energy range the Higgs self-coupling constant 
running equation is 



levr" 



,d\ 
'di 



- ( Igf + 9gi ) A + 



9/3 



gt + ^g!gl + gt]+4YsX-4H 



4 V25 



(40) 
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with 

= lY {3 (rjr,)' + 3 (rjr,)' + {y^y^Y] . (4i) 

The mass of the Standard Model Higgs boson is given by M|- = X{4>ws)'^ aiid, for 
definiteness, we take Mh — 115 GeV at weak scale. 

In order to run the renormalisation group equations down to 1 GeV, we use the 
following initial values: 

U{1) : giiMz) = 0.462 , (?i(Mpi,nck) = 0.614 , (42) 

SU{2) : g2iMz) = 0.651 , (?2(Mpi,nck) = 0.504 , (43) 

SU{3) : gs{Mz) = 1.22 , 53(Mpianck) = 0.491 . (44) 



5 Method of numerical calculation 

According to our philosophy - at the Planck scale, all coupling constants are complex 
numbers of order unity - we evaluate the product of mass-protecting Higgs VEVs required 
for each mass matrix element and provide it with a random complex number of order one 
as a factor. This means that we assume essentially maximal CP violation in all sectors, 
including the neutrino sector. Since the exact values of all these coupling constants are 
not known, our model is only able to provide its results order of magnitudewise. 

In this way, we simulate a long chain of fundamental Yukawa couplings and propagators 
making the transition corresponding to an effective Yukawa coupling in the Standard 
Model. In the numerical computation we then calculate the masses and mixing angles 
time after time, using different sets of random numbers and, in the end, we take the 
logarithmic average of the calculated quantities according to the following formula: 

M=exp(E^) • (45) 

Here (m) is what we take to be the prediction for one of the masses or mixing angles, rrii is 
the result of the calculation done with one set of random number combinations and A'^ is 
the total number of random number combinations used. The calculations are done using 
a Monte Carlo method by putting in for the order of unity couplings random complex 
numbers of order unity. Strictly speaking, we just put such random number factors on the 
different mass matrix elements. We then interpret the matrix elements estimated as the 
products of a random number of order unity and the Higgs VEVs (suppression factors) 
as the running Yukawa couplings at the scale where the Higgs fields W, T, p, uj and x 
roughly all have their VEVs namely about one order of magnitude below the Planck scale. 
We ignore running of couplings over this single order of magnitude as not important. We 
then run down case after case of random numbers the couplings to get the observable 
quantities at the experimental scale, which we take to be 1 GeV as convention so that our 
masses are running ones at = 1 GeV. At first we get the different results for each set of 
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random numbers, but then we average over the logarithm of these quantities and present 
the exponents of the average log's, as "logarithmic averages". In estimating the quality 
of our fits of the adjustable VEVs we define a quantity which we call the goodness of fit 
(g.o.f.) reminiscent of the chi-square using again the logarithms 

2 



g.o.f. = J2 



In 



(m) 



m, 



exp ^ 



(46) 



It should be kept in mind that, in our calculations, we use the renormalisation group 
/^-functions to run Yukawa couplings down to the experimentally observable scale 1 GeV. 
This is because we took the charged fermion masses to be compared to "measurements" 
at the conventional scale of 1 GeV, except for the top quark. We used the top quark pole 
mass instead ISIl: 



/ 4a,(M)' 
where we set M = 180 GeV as an input, for simplicity. 



(47) 



5.1 Mass matrices 



With the system of quantum numbers in Table |2.3| one can easily evaluate, for a given 
mass matrix element, the numbers of Higgs field VEVs of the different types needed 
to perform the transition between the corresponding left- and right-handed Weyl fields. 
The results of calculating the products of Higgs fields needed, and thereby the order of 
magnitudes of the mass matrix elements in our model, are presented in the following mass 
matrices: 



the up-type quarks: 



V2 



to" 



the down-type quarks: 



s) 



V2 



1 lytyt 



iy2(Tt)4 




the charged leptons: 



V2 

the Dirac neutrinos: 



M, 



D 



V2 



cu6(pt)3|y(Tt)2 
\u\p^)\W^fT^ 



(pt)3|ytT2 

1^ (pt)3H^tTt;^t 



(cut)VW^(T^)' 
iy(Tt)2 

(oj^fp^W^T^ 



t^3„3 



WT 



(o;t)VW^tT2x~ 
lytTt 



(4J 



(49) 



(50) 



(51) 
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and the Majorana (right-handed) neutrinos: 



R 



(pt)3;^t/2 ^1 1 

V (P^)V2 1 X 



(52) 



In order to get the true model matrix elements, one must imagine that each matrix 
element is provided with an order of unity factor, which is unknown within our system 
of assumptions and which, as described above, is taken in our calculation as a complex 
random number, later to be logarithmically averaged over as in eq. (^51) . 

The off-diagonal elements of the right-handed neutrino mass matrix (eq. ^2|) are di- 
vided by a factor 2 because the symmetric (Majorana) mass matrix gives rise to the same 
off-diagonal term twice, i.e., we avoid the overcounting of the corresponding Feynmann 
diagrams. However, the element which couples with only one Higgs field should not be 
multiplied by an extra factor 1/2. 

Note that the quantum numbers of our 6 Higgs fields are not totally independent. In 
fact there is a linear relation between the quantum numbers of the three Higgs fields W, 
T and x- 

Qy, = 3Qw-9QT , (53) 

where the 6 components of the charge vector Q correspond to the 6 columns of Table |2.3| . 
Thus the Higgs field combinations needed for a given transition are not unique, and the 
largest contribution has to be selected for each matrix element in the above mass matrices. 

Furthermore, we should mention that we do not assume Supersymmetry, so that we 
can make use of an expectation value of a Higgs field and of the Hermitian conjugate one 
without any restriction; they are numerically equal. If we supersymmetrised the model 
we would need to double the Higgs fields and we would loose predictive power, because 
we would get more parameters. 



5.2 Renormalisation group effect on Dirac neutrino Yukawa cou- 
pling 

In the previous work |1|] we ignored the running of the Yukawa couplings for the neutrino 
sector from the Planck scale down to the see-saw scale and only used the running of 
the effective dimension five operator (eq. ^9|), but in the present article we use also the 
running of the Dirac Yukawa couplings for the neutrinos above the see-saw scale according 
to eq. (|35|). By putting in the various values we may estimate that the most important 
term is the top- Yukawa containing term, Ys, which makes the overall size of the neutrino 
Dirac Yukawa matrix increase towards the ultraviolet. Since our predictions are a priori 
made at the Planck scale, this means that inclusion of this running makes the predictions 
for these Yukawa couplings a bit - about a factor a/2 or less - smaller at lower energy 
scale. 
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6 Results for the quantities of quarks and leptons 



Using the three charged quark-lepton mass matrices and the effective neutrino mass 
matrix together with the renormahsation group equations we made a fit to all the fermion 
quantities in Table || varying just 6 Higgs fields VEVs. We averaged over = 10, 000 
complex order unity random number combinations (see eq. HSI). These complex numbers 



are chosen to be the exponential of a number picked from a Gaussian distribution, with 
mean value zero and standard deviation one, multiplied by a random phase factor that 
has smoothly distributed phase. We varied the 6 free parameters and found the best 
fit, corresponding to the lowest value for the quantity g.o.f. defined in eq. (^6|), with the 
following values for the VEVs: 



■'WS 



) = 246 GeV , = 1.23 x 10'" GeV , {lu) = 0.245 



(p) = 0.256 , (W) = 0.143 , (T) = 0.0742 , (x) = 0.0408 , (54) 

where, except for the Weinberg-Salam Higgs field and {4>b-l), the VEVs are expressed in 
Planck units. Hereby we have considered that the Weinberg-Salam Higgs field VEV is 



already fixed by the Fermi constant. The results of the best fit, with the VEVs in eq. (54), 
are shown in Table ^ and the fit has g.o.f. = 3.38 (see the definition in eq. ^). 

We have 11 = 17 — 6 degrees of freedom - predictions - leaving each of them with a 



logarithmic error of y 3.38/11 ~ 0.55, i.e., we can fit all quantities with a typical error of 

a factor exp (^^3.38/11^ = 1.74 of the experimental value. 

Unlike in older versions of the model, the first and second family sub-matrix of M^, is 
now dominantly diagonal. In previous versions of the model this submatrix satisfied the 
order of magnitude factorisation condition (M£,)i2 • (Mo)2i ~ (M£,)ii ■ (M£,)22; thus the 
down quark mass rrid received two contributions (off-diagonal as well as diagonal) of the 
same order of magnitude as the up quark mass m^. This extra off-diagonal contribution 
to rrid of course improved the goodness of the fit to the masses of the first family, since 
phenomenologically rrid ~ 2 m„. However, in the present version of the model with the 
uj and p Higgs fields, the off-diagonal element {Mj^)2i becomes smaller and we are left 
with a full order of magnitude degeneracy of the first family masses, even including the 
down quark. Our best fit values for the charm quark mass and the Cabibbo angle Vus 
are smaller than in our previous fits to the charged fermion masses |TT|. Nonetheless, 
as mentioned above, our present best fit agrees with the experimental data within the 
theoretically expected uncertainty of about 64% and is, therefore, as well as can expect 
from an order of magnitude fit. 



6.1 Neutrino oscillation parameters 

The charge current interactions from the Sudbury Neutrino Observatory (SNO) ||12| have 
provided an important signal confirming the existence of the solar neutrino mass [p!3| , |1^, 
|15|, |l^, |l^ : SNO detected a flux of z/^ neutrino or Ur neutrino among solar neutrinos after 
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traveling from the core of the Sun to the Earth. Combination of the SNO results with 
previous measurements from other experiments confirms the standard solar model ||T8| , 
whose predictions of the total flux of active ®B neutrinos in the Sun agree with the 
SNO and Super-Kamiokande [17 data. Furthermore, the measurement of the and 
hep solar neutrino fluxes shows no significant energy dependence of the electron neutrino 
survival probability in the Super-Kamiokande and SNO energy ranges. In fact, global 
analyses p9| , pO| , |2l| , p^] of solar neutrino data, including the first SNO results and the 



day-night effect pSf, which disfavoured the SMA-MSW solution at the 95% C.L., have 
confirmed that the LMA-MSW solution gives the best fit to the data and that the SMA- 
MSW solution is very strongly disfavoured and only accepted at the 3a level. 



Furthermore, the combination of the results from atmospheric neutrino experiments p4 



and the CHOOZ reactor experiment constrains the first- and third- generation mix- 
ing angle to be small, i.e. the 3a upper bound is given by tan^ ^chooz ~ 0.06. This limit 
was obtained from a three flavour neutrino analysis (in the five dimensional parameter 
space - ^0, ^chooz, Oatm, Arftg and Aml^^), using all the solar and atmospheric neutrino 
data and based on the assumption that neutrino masses have a hierarchical structure, i.e. 





Fitted 


Experimental 


rriu 


5.2 MeV 


4 MeV 


rrid 


5.0 MeV 


9 MeV 


me 


1.1 MeV 


0.5 MeV 


rric 


0.70 GeV 


1.4 GeV 


rus 


340 MeV 


200 MeV 


mf. 


81 MeV 


105 MeV 


Mt 


208 GeV 


180 GeV 


rrib 


7.4 GeV 


6.3 GeV 




1.11 GeV 


1.78 GeV 


^us 


0.10 


0.22 


Veb 


0.024 


0.041 


Vub 


0.0025 


0.0035 


Am| 


9.0 X 10-^ eV^ 


4.5 X 10"^ eV^ 




1.8 X 10-3 eV^ 


3.0 X 10-3 eV^ 


tan^ 9q 


0.23 


0.35 




0.83 


1.0 


tan ^chooz 


3.3 X 10-2 


< 2.6 X 10-2 


g.o.f. 


3.38 





Table 2: Best fit to conventional experimental data. All masses are running masses at 
1 GeV except the top quark mass which is the pole mass. Note that we use the square 
roots of the neutrino data in this Table, as the fitted neutrino mass and mixing parameters 
(m), in our goodness of fit (g.o.f.) definition, eq. (|4^). 
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In Table |] presented solar neutrino data, however, come from a global two flavour 
analysis, which means that the first- and third-generation mixing angle is essentially put 
equal to zero, i.e., the dependence of ^13 on the solar neutrino parameters have been 
ignored. In principle, of course, we have to fit to neutrino parameters from a three flavour 
analysis. 



The global three flavour analyses have been done by several authors |^B|, ^ ^ ^ and 



they showed a significant influence of the non-zero CHOOZ angle on the solar neutrino 
mass squared difference and mixing angle. In |^ the relatively large solar neutrino mass 
squared difference lying in the LMA-MSW region (with the condition AmQ>,10~^ eV^), 
the solar mixing angle and the CHOOZ reactor experiment data were analysed using the 
three flavour analysis method. Their works tell us that if the CHOOZ angle becomes big- 
ger than zero then the solar mixing angle becomes smaller. This effect is more significant 
for the larger Attt-q values. Because of this correlation, our fit to the neutrino data should 
be somewhat better than that suggested by the g.o.f. value; even including the CHOOZ 
angle our neutrino fit is extremely good. 

Experimental results on the values of neutrino mixing angles are usually presented in 
terms of the function sin^ 26 rather than tan^ 6 (which, contrary to sin^ 26, does not have 
a maximum at ^ = 7r/4 and thus still varies in this region). Transforming from tan^ ^ 
variables to sin^ 26 variables, our predictions for the neutrino mixing angles become: 

sm'^26Q = 0.61 , (55) 
sin2 2^,tm = 0.99 , (56) 
sin2 2^ehooz = 0.12 . (57) 

We also give here our predicted hierarchical neutrino mass spectrum: 

mi = 9.8 X 10~^ eV , (58) 
m2 = 9.6 X 10"^ eV , (59) 
7713 = 4.4 X 10~^ eV . (60) 



Compared to the experimental data these predictions are excellent: all of our order 
of magnitude neutrino predictions lie inside the 95% C.L. border determined from phe- 
nomenological fits to the neutrino data, even including the CHOOZ upper bound. On the 
other hand, our prediction of the solar mass squared difference is a factor of 2 larger than 
the global fit data even though the prediction is inside of the LMA-MSW region, giving 
a contribution to our goodness of fit ~ 0.12. 

Our CHOOZ angle also turns out to be about a factor of ^/2 larger than the exper- 
imental limit at 90% C.L., corresponding to another contribution of g.o.f. ~ 0.014. In 
summary our predictions for the neutrino sector agree extremely well with the data, giv- 
ing a contribution of only 0.25 to g.o.f. while the charged fermion sector contributes 3.13 
to g.o.f. 
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Note that the value of {4>b-l) presented in Table ^ does not fit very well the atmo- 
spheric neutrino mass squared difference, being about factor a/2 less than the best fit 
reported by Super-Kamiokande collaboration p4|. That means that our model predicts 
"relative" degenerated mass spectra in the neutrino sector. However, if we force it to fit 
the value of mass squared difference as Aml^^^ = 2.5 x 10^'^ eV^ and just arrange the 
scale, {(Pb-l), of the mass squared differences - it will not significantly influence the ratio 
of mass squared differences - thus we get AtjIq = 1.2 x 10~^ eV which lays still in the 
allowed region of the global fit of the solar neutrino analysis: 



■>B-L 



Am 



atm 



corrected 
corrected 
corrected 



1.0 X 10^° GeV , 
2.5 X 10"^ eV^ , 
1.2 X 10"^ eV^ . 



(61) 
(62) 
(63) 



and the CHOOZ angle does not change tan ^^chooz = 3.3 x 10 in this small deviation of 
the see- saw scale. 



6.2 CP violation 

Since we have taken our random couplings to be - whenever allowed - complex we have 
order of unity or essentially maximal CP violation so that a unitary triangle with angles 
of order one is a success of our model. After our fitting of masses and of mixings we can 
simply predict order of magnitudewise the CP violation in e.g. K^—K^ decay or in CKM 
and MNS mixing matrices in general. 

The Jarlskog area Jcp provides a measure of the amount oi CP violation in the quark 



sector ||3^ and, in the approximation of setting cosines of mixing angles to unity, is just 



twice the of the unitarity triangle: 

Jcp = Vus Vcb Vub sin 5 , (64) 

where 5 is the CP violation phase in the CKM matrix. In our model the quark mass 
matrix elements have random phases, so we expect 5 (and also the three angles a, (3 
and 7 of the unitarity triangle) to be of order unity and, taking an average value of 
I sin (5 1 ~ 1/2, the area of the unitarity triangle becomes 

Jcp ~ ^ Vus Vcb Vub ■ (65) 

Using the best fit values for the CKM elements from Table |^, we predict Jcp ~ 3.0x 10~^ to 
be compared with the experimental value (2 — 3.5) x 10~^. Since our result for the Jarlskog 
area is the product of four quantities, we do not expect the usual ±64% logarithmic 
uncertainty but rather ±a/4-64% = 128% logarithmic uncertainty. This means our result 
deviates from the experimental value by ln( gQ^|^[|_6 )/1.28 = 1.7 "standard deviations". 

The Jarlskog area has been calculated from the best fit parameters in Table |^, it is 
also possible to calculate them directly while making the fit. So we have calculated Jcp 
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for = 10, 000 complex order unity random number combinations. Then we took the 
logarithmic average of these 10, 000 samples of Jcp and obtained the following result: 

= 3.0 X 10^6 , (66) 

in good agreement with the values given above. 



6.3 Neutrinoless double beta decay 

Another prediction, which can also be made from this model, is the electron "effective 
Majorana mass" - the parameter in neutrinoless beta decay - defined by: 



3 

^ei 

i=l 



(67) 



where rrii are the masses of the neutrinos Ui and Uei are the MNS mixing matrix elements 
for the electron flavour to the mass eigenstates i. We can substitute values for the neutrino 
masses from eqs. ( p8| - p0| ) and for the fitted neutrino mixing angles from Table ^ into 
the left hand side of eq. (^). As already mentioned, the CP violating phases in the MNS 
mixing matrix are essentially random in our model. So we combine the three terms in 
eq. (^) by taking the square root of the sum of the modulus squared of each term, which 
gives our prediction: 

|(m)| ^ 2.5 X 10"^ eV . (68) 

In the same way as being calculated the Jarlskog area we can compute using A^ = 
10, 000 complex order unity random number combinations to get the | {itl)\. Then we took 
the logarithmic average of these 10,000 samples of \{m)\ as usual: 

|(m)| = 3.4 X 10-3 eV . (69) 

This result does not agree with the central value of recent result - "evidence" - from the 
Heidelberg- Moscow collaboration [pl[]. 



7 Baryogenesis via lepton number violation 

Now we have a good model which predicts orders of magnitude for all the Yukawa 
couplings including the see-saw particles, so it is natural to ask ourselves whether this 
model can predict also the right amount of Yb - the ratio of the baryon numbers density 
relative to the entropy density - using the Fukugita-Yanagida mechanism . According 
to this mechanism the decays of the right-handed neutrinos by CP violating couplings 
lead to an excess of the B — L charge (baryon number minus lepton number) , the relative 
excess in the decay from Majorana neutrino generation number i being called e^. This 
excess is then immediately - and continuously back and forth - being converted partially 
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to a baryon number excess, although it starts out as being a lepton number L asymmetry, 
since the right-handed neutrinos decay into leptons and Weinberg-Salam Higgs particles. 

It is a complicated discussion to estimate to what extend the B — L asymmetry is 
washed-out later in the cosmological development, but in our approximation below we 
agree within a factor 3 with the baryon number excess left to fit the Big Bang devel- 
opment at the stage of formation of the light elements primordially (nuclearsynthesis). 
The "experimental" data of the ratio of baryon number density to the entropy density is 
obtained by recent measurement of cosmic microwave background radiation p3| : 



B 





' +1.5\ 


exp ( 









X 10 



-11 



(70) 



Recently, baryogenesis calculations had been done by several authors 
ent models based on this scenario, and they used "usual" range of Yb'- 



Y, 



B 



:i.7-8.1) X 10"^^ . 



exp 



using differ- 
(71) 



In the following subsection, we review briefly the already known Fukugita-Yanagida 
mechanism. 



7.1 Fukugita-Yanagida scenario for lepton number production 

The SU{2) instantons ||3^, rather say, Sphaleron |^ guaranteed the rapid exchanges 



of the minus baryon number and lepton number in which though B — L is conserved in 
the time of Big Bang, even when the temperature was above the weak scale. In fact 
the three right-handed neutrinos decay in the B — L violating way, at the temperature 
about the see-saw mass scales. This means that the baryon number is violated. From 
our assumption - all fundamental coupling constants are of order of unity and we treat 
them as complex random numbers at the Planck scale - it is clear that the model has not 
only C violation but also CP violation. Out-of-equilibrium condition comes about during 
the Hubble expansion due to the excess of the three type of the right-handed Majorana 
neutrinos caused by their masses. 

These statements lead to that all Sakharov conditions are fulfilled in our model: (1) 
baryon number violation, (2) C and CP violation, (3) departure from thermal equilibrium. 



7.2 Entropy density in cosmology 

In order to investigate the quantity, Yb, we need the expression for the entropy of Planck 
radiation which is given by 

., = ^^T^ (72) 

where g^,i is the total number of effectively massless degrees of freedom of the plasma 
at the temperature of the right-handed neutrino, and the index i denotes the number of 
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copiously existing right-handed neutrinos at the time in question. For the estimate of 
the excess of baryon number coming from a certain right-handed neutrino i, we should 
compare to the entropy density at the temperature being approximately equal to the mass 
of the corresponding right-handed neutrino. 

The are obtained as follows: There are 14 bosons and 45 well-known Weyl fermions 
plus i Majorana particles when the temperature is high enough compared to 10^ GeV: 

j=bosons j=fermions 

7 7 
= 28 + - ■ 90 + --i 

8 4^ 

Standrad Model see-saw particles 



108.5 
110.25 
112 



i = 1 

t = 2 , (73) 
i = 3 



here Tj denotes the effective temperature of any species j. When we have coupling as at 
the stage discussed between all the particles Tj = T. 



7.3 CP violation in decays of the Majorana neutrinos 

A right-handed neutrino, N^i, decays into a Weinberg-Salam Higgs particle and a left- 
handed lepton or into the CP conjugate channel: These two channels have different lepton 
numbers ±1. If there were no CP violation they would have the same branching ratio. 
However, all Yukawa couplings are complex order unity random numbers, therefore the 
partial widths do not have to be equal in the next-to-leading-order of perturbation theory: 
CP violation in the decay of right-handed neutrinos. 

Defining the measure for the CP violation 



- - - ■ (74) 

' Ea,/3r(iv^,^£"0^s) + Ea,/3r(iv,,^£"0^^,) ' 

where F are N^i decay rates (in the rest frame), summed over the neutral and charged 
leptons (and Weinberg-Salam Higgs fields) which appear as final states in the N^i decays 
one sees that the excess of leptons over anti- leptons produced in the decay of one iV^j is 
just ej. Now we shall calculate in perturbation theory: The total decay rate at the tree 
level (Fig. ||(a)) is given by 

((M^)tMj).. 

Ttv, = Fjv^^ + F^^£- = — -^^y- — Mi , (75) 

where can be expressed through the unitary matrix diagonalising the right-handed 
neutrino mass matrix Vr: 

= Vr , (76) 
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(a) 



-'ws 



^ws 



■'ws 



Figure 2: Tree level (a), self-energy (6) and vertex (c) diagrams contributing to heavy 
Majorana neutrino decays. 



V^MRMiVR = diag{MlMlMi] 



(77) 



The CP violation in the Majorana neutrino decays, e^, arises when the effects of loop are 
taken into account, and at the one-loop level, the CP asymmetry comes both from the 
wave function renormalisation (Fig. 0(&)) and from the vertex (Fig. 0(c)) I^S], |40[] : 




(78) 



where the function, f{x), comes from the one-loop vertex contribution and the other 
function, g{x), comes from the self-energy contribution. These e's can be calculated in 
perturbation theory only for differences between Majorana neutrino masses which are suf- 
ficiently large compare to their decay widths, i.e., the mass splittings satisfy the condition, 
|M - MA > IF. - FJ: 



fix) 
9{x) 



1 - (1 + x) In 



1+x 



X 



1 — X 



(79) 
(80) 



The function f{x) has for large x "tricky" cancellations, so for numerical calculation we 
use following approximation in the region x ^ 1, i.e., hierarchical case in the right-handed 
sector: 

fix) ^ . (81) 



2^ 



7.4 Off-diagonal dominant matrix elements in the see-saw neu- 
trino mass matrix 

The main modification in the model - and not only in the calculation - compared to the 
latest version is that we changed the quantum number assignment of the see-saw scale 
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fitting Higgs field, (f)B-L, so that it instead of giving directly and thus unsuppressed the 
(3, 3) mass matrix element, i.e., a diagonal one, it gives in the present version unsuppressed 
size to the (2, 3) and (3, 2) mass matrix elements in this see-saw neutrino mass matrix 
(eq. ^). By having such off-diagonal dominance one expects a large mixing of the second 
and third generation at least in the see-saw mass matrix and expect to get easily a large 
CHOOZ angle tan^ 6'chooz- However, as we have already presented in the section]^, we still 
do fit the CHOOZ angle limit satisfactorily. 

The gain by having this off-diagonal dominance is, however, crucial for the baryon 
number production in the early Universe via the Fukugita-Yanagida mechanism: By hav- 
ing now the two heaviest see-saw neutrinos almost degenerate in mass the interaction 
between these two levels become much enhanced and this causes the CP violation as 
expressed by the e^'s for the two heaviest Majorana neutrinos to be bigger, as is easily 
understood by noticing the factor f{x) + g{x) occurring in the expression (eq. ^), which 
is in first approximation a suppression factor given by the ratio of the see-saw neutrino 
masses . With the mass difference becoming small compared to the masses them- 

selves the term g{x) can even become big compared to unity, so that one might see the 
possibility of these normally suppressing factors causing an enhancement. However, there 
are in the case with two such close Majorana neutrinos a strong cancellation between the 
amounts of B — L produced in excess by these almost mass degenerate particles decay- 
ing. Nevertheless we can avoid significantly the suppression by having such degeneracy, 
therefore the baryon number production increases. 



8 An approximation to the wash-out effect 

The excess B — L produced at first gets partly destroyed by the same or other B — L 
violating processes. There are several processes/Feynmann diagrams that can lead to 
B — L violation and thus wash-out an excess. Most important at the temperature scale of 
the right-handed neutrino masses are the resonance scattering of a Weinberg-Salam Higgs 
towards a left-handed lepton producing one of the Majorana neutrinos as a resonance. 
The rate of forming resonances this way is proportional to the width Fj of the see-saw 
neutrino functioning as resonance from detailed balance. Since the time scale involved is 
given by the Hubble expansion at the time of the temperature being equal to the mass of 
the see-saw neutrino, a crucial parameter for the wash-out effect via the resonance process 
is 

K,^I^ = ^"-r ((Ag)'Ag). (,^i_2_3)_ ,32) 

2H T.M. i.66{^,„)^8i9,'f Af, '■I 

where F, is the width of the flavour i Majorana neutrino, Mj is its mass and g^^i is the 
number of degrees of freedom at the temperature Mj. 

As the time goes in early cosmology, the heaviest see-saw neutrino goes out-of-thermal- 
equilibrium first and deploys its excess of B — L, then the next and so on. Thus in the 
hierarchical case - but, our model is not this case in as far as the two heaviest see-saw 



24 



neutrinos are approximately degenerate in mass - the excess produced by the hghtest see- 
saw neutrino will roughly only be washed out by itself but not by the heavier neutrinos, 
since the latter resonances are not reachable at the temperature at that time. However, 
the excess from the decay of the heavier see-saw neutrinos can be washed out by the 
lighter one(s), and if some are degenerate they may wash out the products of each other. 
Now the flavours of the left-handed leptons produced in excess in the decays are not the 
same for the three different see-saw neutrinos. They are, however, also not orthogonal 
states because they are mixing almost maximally. This means that wash-out due to one 
see-saw neutrino of the excess caused by another one is suppressed compared to what the 
wash-out rate would be when a see-saw neutrino washes out its own production of excess. 
We shall take into account such suppressions by constructing some effective Ki called K^g^ 
which should mean that value of Ki to use as if we had the excess of the z/^j washed out 
by itself, so as to obtain in practice the effect of the wash-out due to the other see-saw 
neutrinos included. The decay products of a given see-saw neutrino are found from the 
Yukawa couplings of this right-handed neutrino to a Weinberg-Salam Higgs and a left- 
handed lepton, and these couplings are proportional to the M^-mass matrix elements. 
We - crudely - shall estimate K^ni expressions in terms of KiS which are the parameters 
for the wash-out of the products of the right-handed see-saw neutrinos itself, where i = 1 
is the lightest and three {i = 3) the heaviest. 

In order to estimate the effective K factors we first introduce some normalized state 
vectors for the decay products: 



1 



ki 



( [M^{Mi)]^. , [M^{Mi)]^^ , [M^{Mi)]J , (83) 



and further define their overlap: 

Then we may take an approximation for the effective K factors: 

i^effl = i^l(Ml) , 

K^,, = K^iM^) + C23 K,{M,) + C21 i^i(Mi) , 
K^,, = Ks{Ms) + C32 K^iM^) + C31 ^i(Mi) . 



(84) 

(85) 
(86) 
(87) 



Strictly speaking, one shall namely not imagine the excesses to flow around as left-handed 
leptons only, because there is with high rate the weak instanton or Sphaleron activity 
converting e.g. lepton number into minus baryon number. Since, however, the instanton 
process at the same time convert one particle from every generation, there remains the 
same excess of one generation over the other even after the instanton-process. Thus we 
expect that the assumed suppressed wash-out eflcct from a different see-saw neutrino 
works anyway. We have for simplicity ignored the scattering processes due to exchange 
processes as candidates for wash-out effects, expecting that at the temperature scales 
relevant they will effectively be small due to cross-section going as the fourth power of 
Dirac neutrino Yukawa couplings rather than the square as does the resonance cross- 
section averaged over energy. 
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8.1 The dilution factors 



When the K^f^^s are not small (not less than one), we have to take into account the dilution 
effects - wash-out effects - for the calculation. Therefore, we define the suppression factor 
Hi to be fraction oi B — L excess produced by the right-handed neutrino number i which 
survives. That is to say, we define so that the resulting relative to entropy density, Sj, 
baryon number density minus lepton number, i.e., Yb-l becomes 



Yb-l = 



i=i 9* 

A good approximation for Kj, the dilution factor, is inferred from refs. [^T], 



10^ 

1 < K < 1 n 
r\< K . < 1 



Ki 



(0.1i^',ffi)2 exp 
0.3 



-^(0.1ir.«,)i 



1 



3 ; 
|5 



2 K^ffi 



(89) 
(90) 

(91) 
(92) 



Note that these dilution factors are not smooth; therefore we used in numerical calcu- 
lations instead of the eq. (|9lD and eq. (|92D the following interpolating redefined dilution 
factor in the range <, K^g^ ^10: 



0^ -f^effi ^$/10 : 



(93) 



8.2 B — L to baryon number conversion 



We have presented all quantities to construct Ib-l in foregoing sections, moreover, 
we should note here that due to the electroweak Sphaleron effect, the baryon number 
asymmetry Yb is related to the baryon number minus lepton number asymmetry Yb-l 0] : 

8A^/ + ANh 



Yj 



B 



22Nf + 13Nh 



Y, 



B-L 



(94) 



where Nf is the number of generations and Nh the number of Higgs doublets. It turns 
out that the final form of the baryogenesis using eq. (^8]) is: 



Yn 



28 
79 



3 
i=l 



9* i 



(95) 
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9 Result of baryogenesis 



The numerical calculation of baryogenesis will be presented in this section. All the 
Yukawa couplings - VEVs of seven different Higgs fields - are obtained by fitting the 
fermion quantities (see section 6). In order to get baryogenesis in Fukugita-Yanagida 
scheme, we have to calculate the following informations, i.e., the see-saw neutrino masses, 
K factors and CP violation parameters (e^'s). At first, we present the see-saw neutrino 
masses then -ft'oHj's using = 10,000 random number combinations and logarithmic 
average: 

Ml = 2.1 X 10^ GeV , (96) 
Ma = 8.8 X 10^ GeV , (97) 
M3 = 9.9 X 10^ GeV , (98) 

and 

(99) 
(100) 
(101) 

The numerical results of our best fitting VEVs also gives the e^'s using same method: 

|ei| = 4.62 X 10"^2 , (102) 
leal = 4.00 x 10"^ , (103) 
lesi = 3.27 X 10"^ . (104) 

The sign of is unpredictable due to the complex random number coefficients in mass 
matrices, therefore we are not in the position to say the sign of e's. It turns out that 
we should calculate directly Yb instead of putting the results from eqs. (|96| ) - (|98|) and 
eqs. ( |102| )-( p!04D using the dilution factor formulae into eq. (p5|). Using the complex order 





= 31.6 


A'cff2 


= 116.2 


-^cff3 


= 114.7 



unity random numbers being given by a Gaussian distribution we get after logarithmic 
average 

= 2.59+225X10"'' ' (105) 



hereby we estimate the uncertainty in the natural exponent according ref. [|rT| to be 
64 % ■ v^TO ~ 200 %. Dominantly the signs of are strongly correlated, therefore, 
we included in the numerical calculation using eq. (^) the correlation of signs effects 
correctly. 

The other remark of the scale dependence for our prediction is on the baryogenesis 
calculation: It is obvious that the quantities, ej, are approximately independent of the 
scale {(pB-b) from the eq. (|7S|). In fact these quantities are influenced by only the ratio 
of the different mass eigenstates of the right-handed neutrino, and in addition through 
the Dirac neutrino Yukawa couplings which are, of course, depended on the see-saw scale 
via renormalisation group equations. However, a small deviation of the see-saw scale can 
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not contribute significantly. The number of degrees of freedom, g^i, does not depend on 
the see-saw scale. However, the dilution factors, rather say, the K^g^ factors are strongly 
influenced by the right-handed scale, i.e., inverse power of the right-handed neutrino 
masses (see eq. ^2]). Strictly speaking, the wash-out effect is proportional to the see-saw 
scale. Using the "right" value of the see-saw scale which we mention in above (eq. |61|) we 
get the corresponding K^g's as following 

(106) 
(107) 
(108) 





corrected 


= 37.60 




corrected 


= 138.8 




corrected 


= 137.1 



The corresponding result of baryon asymmetry is 



Yn 



corrected 



+ 13.3 11 
-1.75 



(109) 



This result is a bit smaller than the result in eq. ( |105| ), however, still in allowed region of 
the "experimental" data in eq. (|7lD. 



10 Proton decay 

It should be mentioned that our non-supersymmetric model passes without any trou- 
ble the test of not predicting proton decay that should have been observed in present 
experiment. Indeed we have in our model desert, except for the see-saw neutrinos and 
associated Higgs, and gauge U{1)b~l fields, up to about one order of magnitude 

below the Planck scale. 

A reasonable suggestion is that the baryon decay causing bosons would have masses 
of the order of the scale around the VEVs of our Higgs fields, which typically are about 
a factor 10 below the Planck scale (see eq. (0)). We should therefore use such a boson 
with mass about 1.2 x 10^^ GeV. That is to say we could essentially use the formula for 
the life time of the proton Tp in simple SU (5) GUT with replaced mass of X-particle by 
1.2 X 10^^ GeV. We could take for instance the expression for the p tt^ channel 

34 / 0.015 GeV^y / My 



"(^-°--) = r(^-.vrOe+) ='"^° -[ 'a, ) " llO^^j ' ^'''^ 

where an is the hadronic matrix element and My the mass of the GUT gauge boson, 
and insert 1.2 x 10^^ GeV. Then we get crudely our prediction assuming the existence of 
bosons which break baryon number at this scale: 

r(p^^Oe+) ^ 10"^^ years . (Ill) 

However, non of the Higgs fields and gauge bosons, which we truly consider in our 
model would be able to cause baryon decay. In fact we have excluded by our assumption 
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about our gauge group the possibility of a gauge boson that could have made transitions 
between quarks and leptons: we excluded, namely, gauge bosons unifying different, irre- 
ducible representations of the Standard Model, and of course - it is well-known - quarks 
and leptons belong to different irreducible representations. Concerning the possibility 
that our scalar fields (Higgs fields) considered be able to cause baryon decay, it should be 
remarked that since they are made to break our gauge group (eq. ||) but to leave unbro- 
ken the Standard Model gauge group. Therefore, they must contain a singlet under the 
diagonal group in order not to break the diagonal SU{3) - QCD gauge group. Then it 
follows that they must under this diagonal SU{3) have triality t = (mod 1) (see below 
eq. (0)). In turn, that means that these Higgs fields are not able to couple to neither two 
quarks nor to a pair of a quark and a lepton. That makes our scalar fields of no practical 
use for baryon decay since to look for three quarks coming together to form one point 
particle at a scale of 10^^ GeV is too much requirement. 

If we keep to truly postulated particles in our model we would have to wait for proton 
decay till we can get sensitive to it from Planck scale masses inserted for My- Taken this 
attitude of only using "the everything" that can be found at Planck scale will lead to a 
10^ times longer proton life time prediction in our model: 



T(, 



(p^wO e+) 



^ 10^^ years . (112) 

Planck ^ ^ 



In any case there is no chance in foreseeable future to observe proton decay according 
to the suggestions of our model - however presumably - this is a point that could be 
modified without changing our model too drastically. 



11 Conclusions 



We have set up a model for calculating all the mass matrices for quarks and leptons - 
here under especially the effective left-handed neutrino mass matrix which is obtained via 
the well-known see-saw mechanism using a set of three see-saw neutrinos, two of which 
are approximately degenerate in mass. We have fit all the masses and mixing angles 
within the accuracy which we theoretically expect to be ±64% (on a logarithmic scale) 
for the quantities for instance quark and lepton masses which are essentially directly 
given by mass matrix elements. In fact even our worst deviations of our predictions 
form experiment are of these to mass matrices simply related quantities: we predict the 
charm quark mass and the Cabibbo angle a factor two lower than experiment, while we 
obtain the electron and the strange mass a factor two larger. On the other hand, the 
squared quantities of neutrino oscillations presented in Table || do not even deviate a 
factor two. The CP violating parameter - Jarlskog triangle area -, Jcp, and the baryon 
number relative to the entropy ratio, Yg, arise as products and ratios of several mass 
matrix elements and thus have larger expected uncertainties, and turn out indeed to 
agree quite within the expected accuracy both being predicted to the low side of "data" . 
In conclusion we fit all the Standard Model parameters related to the mass matrices within 
the expected order of magnitude accuracy, i.e., we fitted 9 masses of quark and charged 
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leptons, 2 mass squared differences for neutrinos, 5 genuinely measured mixing angles, tlie 
Jarlskog triangle area and the baryon number to entropy density ratio. In addition, we 
managed to avoid violating the bounds for proton decay life time and neutrinoless double 
beta decay. Most interestingly we predicted the CHOOZ mixing angle - the mixing angle 
of the electron neutrino with the heaviest eigenmass neutrino - only about 15% above the 
experimental limit (completely allowed order of magnitudewise). In this way we predicted 
successfully 18 genuinely measured parameters not given by the Standard Model itself, 
and that we did with only 6 adjustable parameters - the vacuum expectation values of 
Higgs fields - in our model. This means that we made 18 — 6 = 12 genuine predictions in 
addition to avoiding limits giving potential problems for most interestingly the CHOOZ 
angle, the proton life time and the neutrinoless double /9-decay, however, we also did not 
predict the rate of yU ^ e + 7 nor similar suppressed processes with any problematic rate. 
If one takes it as we mentioned in the earlier article that the smallness of the fine 
structure constants in a model with our gauge group is understandable as being really of 
order unity except for the An in = gf/An and the corrections from the break down to 
diagonal subgroup, we can claim that we have understood with our model all the orders 
of magnitudes of the couplings and parameters of the Standard Model in terms of only 
7 vacuum expectation values (if we should claim the VEV of Weinberg-Salam Higgs 
as understood we would, of course, have to include it simply as a seventh parameters 
because we have not explained the hierarchy of the weak to Planck scale). 

We did obtain this impressive fit order of magnitudewise in a model with the gauge 
group eq. (P and by a very data inspired choice of the gauge quantum numbers of our seven 
Higgs fields. We have, however, good reasons to claim that the precise detailed structure 
of our gauge group is not important, especially in as far as we in reality only needed 
to use the abelian part of the extension of the Standard Model group. The presumably 
only important thing is that the gauge group is sufficiently large, so that it separates the 
various Weyl components of the quarks and the leptons, in order to accomplish the various 
mass matrix element suppressions needed. On the other hand, the gauge group should 
not be so large as to include SU{5) or other Grand unification gauge group, since that 
tends to give exact mass relations which become a severe problem to avoid, except for the 
T — b mass relation. Our gauge group were the biggest one represented on the supposed 
three heavy neutrinos and the 45 already observed Weyl quark and lepton components 
and avoiding such unwanted unifications. However, such a large gauge group is far from 
being needed although it must be considered the suggestion of the present article that a 
relatively large group mainly abelian or, at least, not too much "unifying" is what could 
fit very well the "small hierarchies" of quark and lepton masses and mixings. 

The choices of the specific quantum numbers for our seven Higgs fields represents an 
opportunity of a sort of discrete fitting in the sense that we adjust these quantum numbers 
to be able to fit the masses and mixings. Historically we developed the quantum number 
proposals making small changes gradually along so as to incorporate new experimental 
data. For example the fields x ci-nd (pB-L were introduced to cope with the neutrino 
oscillations, and previous Higgs fields (called S and ^) were replaced by the fields p and 
(jj in order to fit the Large Mixing Angel MSW solar neutrino solution, which became 
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favoured strongly during the development period of our model. 



Finally the new step of development in the present article is that the quantum num- 
bers of the field (pB-L got modified to make the two heaviest see-saw neutrinos become 
degenerate in mass only deviating by a relatively small mass difference of the order of the 
suppression factor associated with the VEV of the field x being small. Thereby we could 
enhance the CP violation in the decay of the heavier see-saw neutrinos, thus achieving 
a larger baryon number. The latter were needed in as far as the previous version of the 
model, which did not have such a degeneracy, gave too little baryon number excess. 

Although in the present model we get somewhat more wash-out effect we still get closer 
to the from Big Bang fitting estimated baryon number density relative to entropy density, 
and we predict it only a factor three (see eq. |105[ ) below the experimental data, and that 
should be counted as only ln(8.5/2.59)/(vT0 ■ 0.64) = 0.59 "standard deviations"! 

11.1 How does the model function? 

A few comments about how the model functions may be on its place: The diagonal 
elements in the four Dirac mass matrices - the two for quarks, the charged leptons, 
and the "neutrino Dirac mass matrix" - are suppressed to just the same degree in all 
these four matrices. This is so because the quantum number exchanges needed to uphold 
these elements are just those of a conventional Standard Model Weinberg-Salam Higgs 
field, however, taken as the family quantum number for that family that corresponds to 
the number of the diagonal matrix element in question. This gives at first an order of 
magnitude degeneracy of quarks and leptons in the same family at the Planck scale of 
running couplings, a result that is good for simulating SU{b) GUT predictions - only order 
of magnitudewise - and for explaining the crude degeneracy of up-quark and electron mass 
(when extrapolated to Planck scale). However, the Charm quark and top-quark masses 
need to be explained as shooting up due to some special mechanism, and we managed to 
get them dominated by off-diagonal elements. Due to that the first family gets its mass 
matrix elements using p and uj together with the Higgs fields already used to give the 
second family masses one tends to get just the same extra factor of the ~ type 
and for mixing angles/off-diagonal elements cop^. There is as a result in such models 
easily achievable factorisation of the mixing angles for the quarks and for the neutrinos, 
respectively, at the Planck scale: 



This going to first family via the second family just getting an extra or cup^ roughly 
is also behind how in our model we get much less hierarchy among the left-handed neu- 
trinos than for the charged quarks and leptons. In fact the lightness of the right-handed 
neutrinos "of first generation" gives a large propagator largely being compensated by the 
corresponding smallness of the first family Dirac mass matrix elements. So the much less 



chooz 



Vub 




(113) 
(114) 
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hierarchy for the neutrinos have a rather natural explanation. This fact supports the 
see-saw mechanism and thereby the Fukugita-Yanagida scheme. 

However, we should admit that the rather large neutrino mixing angles in our model 
are obtained by making it possible to fit them to be of order unity. We use e.g. for 
the achievement of the large atmospheric neutrino mixing angle the adjustment of the x 
vacuum expectation value so as to be not much different from that of T. The solar mixing 
angle is made order unity by having p ~ with about same vacuum expectation values. 

Further mass relations at the Planck scale are in our model: 

Vus = e^ ^ (^o)-^ (—Y , (115) 

\ms/ 

ml ~ ms mc mt . (116) 



Concerning the baryon number production it should be said that due to the much 
larger CP violation, i.e., €2 and £3, in the decay of the two heavier and almost degenerate 
see-saw neutrinos than in the decay of the lightest one, it is the heavy neutrinos that 
deliver the main/dominant contribution to the surviving {B — L) excess and thus to the 
baryon number. 

This of course means that it is important for our success that there were no inflation 
after the era of the heavy pair of degenerate see-saw neutrinos to attenuate the B — L 
produced and add re-heating entropy diminishing Yb-l (or Yg). It is also important for 
such success that we had already at this era roughly the number of see-saw neutrinos 
corresponding to thermal equilibrium as we used in our calculation. That means that 
inflation eras should already have been recovered from at the time when the temperature 
reached down to 10^° GeV (see eqs. |96|-p8D. 



11.2 Developments in calculational method 

It should be mentioned that we have made the present numerical calculations taking into 
account the renormalisation group running effects of the Yukawa couplings as well as the 
effective dimension five operator corresponding to two Weinberg-Salam Higgses interacting 
with two left-handed neutrinos giving the neutrino masses observed. Compared to our 
previous work we included the running between the Planck scale and the see-saw scale of 
the Dirac neutrino Yukawa couplings. 



11.3 On the accuracy 

At first one might wonder if such runnings are needed when we at the end in principle 
get only order of magnitude predictions. However, they can give easily factors of three (or 
more) and we really have managed to get agreeing results with accuracy corresponding 
to that the simple mass matrix elements are uncertain only by the theoretically expected 
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±64% logarithmically, so indeed it is important to include corrections that could be of 
the order of a factor two. It is quite remarkable that we are able to work/fit with this 
good accuracy in a model that is in principle only order of magnitude. If it is not just 
because we had too many discrete details, i.e., quantum numbers of Higgs field, to adjust 
it should mean that there is indeed in nature couplings that are a priori of order unity in 
a way that in praxis means that one can count on them being of that order to the rather 
high accuracy we achieved. With this our ±64% accuracy in mind we may remark that 
we have in spite of working only "with orders of magnitude" gotten rather close to the 
experimental uncertainties in the quantities which we fit: For the very light quarks up, 
down and strange quarks one can get variations form one lattice computation to the other 
one that can "well" reach our uncertainty and the neutrino parameters are typically also 
uncertain with such uncertainties. 

Of course one can still be more ambitious than we were in the present article and 
hope for understanding the well measured mass ratios of the charged leptons, but there is 
not at all as much yet to fit as if the quarks and neutrons masses and mixings had been 
equally well measured and thus an order of magnitude fitting is at the present time using 
quite a significant bit of the data that can possibly help us to dig behind the Standard 
Model. 

At the end let us admit that our ambitious model does not have any candidate in it for 
dark matter nor for the LSND neutrino anomaly. That would require totally new elements 
in our model, as for instance a hidden sector. Combining our model with Supersymmetry 
would enforce us to double all the Higgs fields except ly-Higgs field, and the model would 
loose much predictive power drastically. 

Note added 

During the completion of this article, the reference appeared which treats also the (3 
function for the Dirac neutrino Yukawa coupling. Our result (eq. ^) does not match with 
their result, with respect to the coefficient of weakhypercharge gauge coupling, because 
of the different notation: we used the GUT notation for gi coupling constant. 
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